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The annihilation channel is implemented into the front form calculations of the positronium spec- 
trum presented in Ref. jlj. The effective Hamiltonian is calculated analytically. Its eigensolutions 
are obtained numerically. A complete separation of the dynamical and instantaneous part of the 
annihilation interaction is observed. We find the remarkable effect that the annihilation channel 
stabilizes the cutoff behavior of the spectrum. 
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I. INTRODUCTION 

Recently, a calculation of the positronium spectrum 
within the front form formalism was presented (!]]. In 
this work, the entire spectrum was calculated, i.e. all 
sectors defined by the kinetic component of the total an- 
gular momentum, J z , were considered. However, the one 
photon annihilation channel was omitted. It is the pur- 
pose of the present paper to include this channel into the 
calculations, which are completed with this step. 

The article is structured as follows. In Section || the 
implications of introducing the annihilation channel into 
the positronium theory are stated. Next, the spectrum 
of the "full" positronium is calculated and the results 
are presented. Section IV deals with the parameter de- 
pendence of the results, i.e. the dependence on cutoff's 
and convergence parameters. A discussion of the results 
follows. The calculation of the matrix elements of the 
effective Hamiltonian is reviewed in Appx. |A|. 



II. THE ANNIHILATION CHANNEL IN 
POSITRONIUM 

It is a topic of its own merit to construct a manageable 
Hamiltonian for the positronium system. One possible 
way was shown in Ref. [Q , where an effective Fock space, 
consisting of two sectors |ee) and \eej), named P- and 
Q-space was considered. One can justify the absence of 
any higher Fock state from the structure of the ap- 
plied formalism of effective interactions B. The general 
formalism is set up for a non-abelian SU(N) gauge the- 
ory. Unlike QED, the one boson state is absent in these 
theories because of color neutrality. Nevertheless, one 
has to take care of the one photon state I7) in QED, and 
in what follows, we demonstrate how this can be done. 

Firstly, it is important to notice not only the differ- 
ences between the QED Table | and the analogous QCD 
table in H, but also the similarities. Some of the graphs 
occuring in QCD are absent in QED. In particular those 
graphs with a three- or four-boson interaction and the 



instantaneous interactions connecting four bosons. But, 
although an additional sector occurs as a first row and a 
first column in the QED table, neither is a change in the 
higher Fock sectors observed, nor is the ordering altered 
in any way. 
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TABLE I. The Hamiltonian matrix for QED. The sectors 
n are numbered starting at zero. The vertex, seagull and 
fork interactions are denoted by V, S, F respectively. Diago- 
nal matrix elements are symbolized by •. 

In addition to the P- and Q-spaces defined above, we 
introduce the iV-space, i.e. the sector containing the |t) 
states, into the system. The corresponding projector is 



all QN 

N:= l(7)n>((7)«l 



(1) 



The whole procedure of subsequent projections of higher 
Fock states onto the remaining Hilbert space of states 
can be carried out like before [jjj until one arrives at a 
(2 x 2) matrix, operating in the N- and P-space rather 
than in the P- and Q-space 
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From Table | one can read off the interaction of the one 
photon state with all other sectors: the vertex interac- 
tion annihilates the photon into an electron-positron pair 
and a fork interaction scatters it into the sector {eej). 
The latter interaction is already contained in the effec- 
tive interaction, Eq. (||), because of the projection of the 
Q-space. 
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FIG. 1. The positronium spectrum including the annihi- 
lation channel in the J z =0 sector. Parameters of the calcula- 
tion: a — 0.3, A = 1.0 m/. The mass squared eigenvalues M 2 
in units of the electron mass rrif are shown as functions of 
the number of integration points N = N\ = A?2- The triplet 
states, especially l 3 Si, are lifted up by the annihilation inter- 
action, whereas the singlet mass eigenvalues stay the same. 
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FIG. 2. The graphs of the positronium model. Effective 
photon lines are labelled by hashed rectangles. The graphs 
of the annihilation interaction are at the bottom line. 

Although we projected the higher Fock sectors onto 



the lower ones to construct the effective interaction up 
to now, one is, of course, free to project the (lower) I7}- 
sector onto the (higher) |ee)-sector, and obtains 



ffeff(w) = Hpp + H 
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(3) 



Of course we do this for convenience; we could just as 
well solve the eigenvalue problem of Eq. (||). The projec- 
tion is depicted in the two graphs in the bottom line of 
Fig. |^. One is the dynamic annihilation graph, the other 
is the corresponding seagull annihilation graph. The lat- 
ter is a P-space graph and was omitted in Ref. El follow- 
ing the gauge principle of TANG B . 



III. SPECTRUM INCLUDING THE 
ANNIHILATION CHANNEL 

The spectrum including the annihilation channel shows 
the expected properties: the singlet eigenvalues remain 
unchanged, while the triplet states change. An essential 
point in the actual calculations is that one has to use the 
same counterterms for the Coulomb singularity as used 
without the annihilation channel. This is due to the fact 
that the one photon annihilation part of the interaction 
has no additional singularity that needs to be taken care 
of numeri call y, because of the simple energy denomina- 
tor, Eq. (A2). We compiled our results in the form of 
binding coefficients in Table 
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FIG. 3. Deviation of corresponding eigenvalues for J z =0 
and J z =l including the annihilation channel as a function 
of the number of integration points for a=0.3, A=1.0m/. 
The graph shows AM 2 ~ M'(Jz = 0) - M 2 (J Z = 1) for the 
states l 3 ^! (A), 2 A P 1 (□), and 2 1 P 1 (a). 
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One notes a slightly larger breaking of rotational sym- 
metry as in the case of missing annihilation channel. The 
triplet ground states of different J z are still degenerate to 
a very good approximation, but the discrepancy is bigger 
than without the annihilation channel. The dependence 
of these discrepancies between corresponding eigenvalues 
for J z =0 and J z =l on the number of integration points 
is shown in Fig. |3[ The behavior of the curves is similar 
to those of the calculations without annihilation chan- 
nel, cf. Fig. 13 of Ref. 0. An additional plot, Fig. |, 
shows the dependence of the rotational symmetry break- 
ing on the cutoff A. The discrepancies of corresponding 
eigenvalues are almost independent of the cutoff A. 
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FIG. 4. Deviation of corresponding eigenvalues for J z —0 
and Jz=l including the annihilation channel as a function 
of the cutoff A for a— 0.3, A=1.0 m/. The graph shows 
AM 2 := Ml{J z =G) - Ml{J z =l) for the states I 3 Si (A), 
2 3 Pi (□), and 2 1 P 1 (o). 



To make a comparison of our results to those of per- 
turbation theory easier, we show in Fig. ^ the values 
for the principal quantum number n=2 for both theo- 
ries graphically. The structure of the two plots is al- 
most the same, only the 2 1 Sq state and the 2 3 Pq state 
are exchanged in our results. This is due to the cut- 
off dependence of the S-state, which is larger than that 
of the other states (cf. Fig. ||(b) and next paragraph). 
We stress that the results of the perturbative calcula- 
tions change considerably, when the next higher order 
in a is considered. For example, the mass squared of 
the triplet 2 3 Si is, according to Fulton and Martin 
@, M 2 (2 3 5i) = 3.9780186070 up to order a 3 Ry , which 
is nearer to our result than the value of perturbative 
calculations up to order 0(a 4 ), displayed in Fig. [Hj. 
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TABLE II. The binding coefficients of the singlet (B B ), 
the triplet (B t ), and the hyperfine coefficient Chf are listed 
for a — 0.3, Ni — 41,^2 = 11. For a comparison, the coef- 
ficients of equal time perturbation theory are also listed, up 
to order C(a 4 ) (ETPT) and to C(a 6 lna) 



IV. PARAMETER DEPENDENCE OF THE 
SPECTRUM 

The convergence of the eigenvalues with growing num- 
ber of integration points N is the same as the case of 
no annihilation channel. To be explicit, the convergence 
of the eigenvalues can be shown to be exponential. Wc 
performed a \ 2 nt to the function 

f{N) = a-bexp{{N-N )/c}. (4) 

and obtain for the singlet ground state 1 ^So a = 
3.9061000 ± 0.0000008, b = (9.6399 ± 0.0193) x 10~ 4 , 
and c = 2.2372 ± 0.0063, with X 2 = 5.7. We did not 
perform the limit N — > oo, because the accuracy of the 
results for N > 20 suffices to compare to other data. 

The cutoff dependence of the positronium spectrum in- 
cluding the annihilation channel is comparable to that of 
the spectrum without it. However, a striking difference 
occurs: the inclusion of the annihilation channel stabi- 
lizes the cutoff dependence of the eigenstates. In partic- 
ular, the triplet ground state in Fig. ||(a) shows only a 
small dependence on the cutoff, when one compares it to 
the behavior of the same state in the calculations without 
annihilation channel |5J]. One can fit these curves with a 
polynomial in log A. The singlet ground state eigenvalues 
are the same as without the annihilation graph and for 
the triplet one obtains 

M t 2 (A) ~ 3.91392 - 0.00029 log A + 0.00015 log 2 A. (5) 

Comparison with Eq. (22) of Ref. Q shows that the de- 
crease of the triplet with log A is suppressed by including 
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the annihilation channel by a factor of 60. Also the ex- 
cited states, n=2, show a different behavior, when com- 
pared to Fig. 13 of Ref. M. Here, only one state shows 
level crossing as A grows large. The eigenvalue of the 
2 3 Pi state, however, depends only weakly on the cutoff. 
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FIG. 5. The spectrum with annihilation channel as a func- 
tion of the cutoff A: ground states (n=l). The parameters 
for the calculation are a = 0.3, J z = 0, Ni = 25, N2 = 21. 
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FIG. 6. The spectrum with annihilation channel as a func- 
tion of the cutoff A: first radial excited states (n=2). The 
parameters for the calculation are a = 0.3, J z = 0, Ni = 25, 
N 2 = 21. 

The values for the binding coefficients B n and the co- 



efficient of the hyperfine splitting Chf are presented in 
Table |l[ The values are correct for a cutoff A ~ 2m f 
when compared to results of perturbation theory up to 
0(a 4 ). However, the effects of higher order correction to 
perturbative calculations are significant for a large cou- 
pling such as a — 0.3. The result of perturbation the- 
ory to C(a 6 lna) for the coefficient Chf is considerably 
smaller than that to order 0(a 4 ). 

Concluding, we state that the cutoff dependence of the 
spectrum is improved as compared with the case of miss- 
ing annihilation channel. This makes it evident that the 
annihilation channel is a necessary part of the theory. 
Note that in muonium, the annihilation state is absent. 
This can lead to problems in interpreting the stabilizing 
feature of the channel ||. 

As a further investigation of the properties of our 
model, we will vary the coupling constant and interpret 
the spectrum. A similar procedure was performed by 
Dykshoorn et AL. [fr||, who studied coupled integral 
equations for QED-bound states in equal-time quantiza- 
tion with a variational ansatz. They calculate masses for 
the lowest eigenstates of positronium with and without 
the annihilation channel and plot them versus the cou- 
pling constant. The prominent feature of their figures is 
the occurrence of a critical coupling at which the masses 
become smaller than zero. We have performed the anal- 
ogous calculations within our approach. It seems at first 
glance [Fig. ||(a)] as if the eigenvalues, after decreasing 
quadratically with the coupling as expected by pertur- 
bation theory, stabilize at a coupling a ~ 1.5. However, 
further investigation shows that this is merely an effect of 
the cutoff dependence of the spectrum. The eigenvalues 
in Fig. ||(a) were calculated with a cutoff A = 1.0 m/, 
which is too small for a coupling constant of a = 1.0, 
since then the Bohr momentum is of the same order as 
the cutoff. Fig. [§](b) shows clearly that there is a crit- 
ical coupling. The masses calculated with a cutoff of 
A = 20 77i f tend to zero at a ~ 0.5. A similar value was 
found in 



V. DISCUSSION 

In the present work we completed in a sense the theory 
of positronium in front form dynamics set up in Ref. [jp . 
The theory was enlarged by including the one photon an- 
nihilation channel. The inclusion is a test for the consis- 
tency of the model. In particular, the conception of an ef- 
fective theory operating with resolvents can be falsified if 
the annihilation channel cannot be implemented without 
special assumptions. Our results show that the imple- 
mentation of the annihilation channel is unproblematic 
in the sense that the same formalism can be applied as in 
the case of the projection of the effective |ee~7)-sector onto 
the electron-positron sector. As an interesting property 
of the annihilation channel we find a strict separation of 
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the instantaneous and the non-instantaneous interaction: 
the seagull interaction is present only in the J z = sec- 
tor, whereas the dynamic graph has non-vanishing ma- 
trix elements only for J z = ±1. Our approach passes 
another test: both graphs yield the same contributions 
to the eigenvalues and consequently rotational symme- 
try is seen also in the spectrum including the annihila- 
tion channel. Moreover, the inclusion of the annihilation 
channel improves the results for the hypcrfinc splitting. 



We stress the point that the implementation of the an- 
nihilation channel completes the investigations of how to 
construct an effective interaction for the electromagnetic 
fermion-antifcrmion system in the meaning of the method 
of iterated resolvents. We have put all effects of higher 
Fock states into an effective |ee7)-sector, as far as the 
spectator interaction is concerned, i. e. the interactions in 
which the photons are not directly involved. The remain- 
der of the interaction relies in the coupling function of the 
vertices. A hint to this conclusion is the logarithmic cut- 
off dependence of our results. It is clear that the coupling 
constant depends on the cutoff and has to be analyzed. 
A future aim, beyond the scope of this work, will be to 
show that the physical results of our model become inde- 
pendent of the cutoff, as soon as renormalization group 
techniques are consistently applied. This is supported by 
the fact that we find a stabilizing effect of the annihila- 
tion channel on the dependence of the spectrum on the 
cutoff: all eigenvalues show slower variation with grow- 
ing cutoff A when the annihilation graph is added, in 
some cases this even prevents level crossings with other 
states. However, in muonium for instance, this channel is 
absent and it remains unclear which interpretation this 
should yield. We therefore conclude that our model is 
correct as long as the vacuum polarization effects are not 
considered. 



The possibility of one boson annihilation is the main 
difference between QED and QCD in effectively trun- 
cated Fock spaces where the dynamic three- or four-gluon 
interaction is not possible but resides in the coupling 
function. The one boson sector exists only in QED be- 
cause of the constraint of color neutrality in QCD. We 
showed that this sector can be projected onto the (al- 
ready) effective electron-positron sector. 
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TABLE III. The positronium spectrum for a = 0.3, 
A = 1.0m/,JVi = N2 = 21. The non-relativistic notation for 
the terms and the binding coefficients B n := (4 — M n )/a 2 
for J z —0 and J z — + 1 including the annihilation channel 
are shown. The discrepancy between the eigenvalues is 
AB n := P„(J z =0) — B n (J z = + 1). The numerical errors arc 
estimated from the difference between the values for maxi- 
mum and next to maximum number of integration points. 
The actual errors are smaller due to the exponential conver- 
gence of the eigenvalues with N. The k numbers in brackets 
are the errors in the last k digits. 
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FIG. 7. Compiled spectra of positronium with different 
J z = —3, —2, . . . , +3 including the annihilation channel. All 
spectra have been calculated with a = 0.3, A = 1.0m/, 
Ni = Na = 21. The mass squared eigenvalues M„ in units 
of the electron mass m^ are shown. The notation for the 
states is 3S+1 l, J J z . The resolution of the plot is inadequate for 
the multiplets. Nevertheless, the numerical degeneracy of the 
three triplet ground states 3 S^ 1 , Si, and 3 5*} becomes very 
clear. 
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FIG. 8. The spectrum as a function of the coupling con- 
stant for A = 1.0 m/. Here, the eigenvalues seem to converge 
to a stable value as a grows large. However, this is just an 
effect of the cutoff dependence of the spectrum: for a larger 
cutoff A = 20 vtij the eigenvalues become negative at a critical 
coupling a c ~ 0.5. 
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FIG. 9. The spectrum as a function of the coupling con- 
stant for A = 20.0 m/. For this large cutoff the eigenvalues 
become negative at a critical coupling a c ~ 0.5. 



APPENDIX A: CALCULATION OF EFFECTIVE 
MATRIX ELEMENTS 

The implementation of the annihilation channel is 
straightforward, but non-trivial. We describe shortly the 
calculation in the (J z =0) and (J z = ± 1) sectors. The 
corresponding diagrams, displayed in the lowest line of 
Fig. |[ vanish in all other sectors because of the helicity 
of the photon: no angular momentum larger than J = 1 
is possible. The functions derived depend on the light- 
cone momenta, (x, k±), and are given in Table M. They 
have to be added to the functions listed in Ref. [ If, Table 
V, to obtain the full Hamiltonian. Following the proce- 
dure described in Ref. (I]] (c/. also Ref. Q), we use for 
the energy denominator the following averaged kinetic 
energy 
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The denominator in the one photon sector is simple be- 
cause the photon has zero mass: 
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(A2) 



Note that this denominator does not depend on the di- 
rections of the vectors k±, k'j_, i.e. on the angles tp, ip' . 

The matrix elements of the vertex interaction V g ^ q q 
(cf. e.g. Ref. |l(i[| ) split up into their three components 
with different helicity factors, are 



l\V A \l) = -m^/P 
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x(l 
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{ee\V B \i) = -^fP 



1 



-k 



(ee\V c \j) = ^—e~ lX2ip x 8 + *l5 



A 2 °+A 3 



(A3) 



The complete matrix elements are obtained by symmetry, 
since (ee\Vi\^) — (^\Vi\ee). To classify the states accord- 
ing to their quantum numbers with respect to the third 
(kinetic) component of the total angular momentum, J 3 , 
we project on sectors of definite J z . It turns out that 



because of the simple energy denominator, Eq. (A2), the 
only dependence on the angles comes from the projec- 
tors on the J z sectors and is proportional to cosjt^ — ip'} 
or sin{ip ~ ip'}. As a result, all matrix elements of the 
dynamic annihilation graph for J z =0 vanish when inte- 
grated over the angles. Only for J z = ± 1 do some of the 
matrix elements survive the integration. 

Because of the combination of matrix elements with 
the projectors onto the J z sectors, two types of functions 
emerge for J z = ± 1: one is independent of the angles, the 
other has a dependence proportional to exp{±2i(<^— <p')}. 
The latter vanishes after angular integration. The helic- 
ity table [V is given to illustrate the helicity dependen- 
cies. It holds for J z = + 1. The analogous table for J 2 = — 1 
is obtained by the operation 

W l3 {J z = + \; A!,A 2 ) = -AiW y (J,=-l; -Ai, — A2). 

The simple kinematics (x e + Xg = 1) of the seagull an- 
nihilation graph result in a constant contribution of this 
graph to the Hamiltonian matrix. It is 



(ee\S\ee) = -2/3^ x <T^ 



C +A' 



(A4) 



Because of its helicity factors, the graph acts only be- 
tween states with 



out : 


in 


Ai.Ai =TT 


A'i, A 2 =U 


K,y 2 =it 


A'i, A 2 =|| 


Ai, A2 


=TT 


Waa 


W AB 


Wac 





Ai, A2 


=U 


Wba 


Wbb 


Wbc 





Ai, A2 


=IT 


Wca 


Wcb 


Wcc 





Ai, A2 


=11 















TABLE IV. Symbolic helicity table for the dynamic an- 
nihilation graph. The functions Wu are identical with the 
expressions Fi listed in Table [v| Here, terms proportional to 
5\j z \fi are omitted. 



The dependence of the effective interaction on the he- 
licities of in- and out-going particles is displayed in the 
form of tables. The following notation is used for func- 
tions of the type F(x,k±;x',k' A ). An asterisk denotes 
the permutation of particle and anti-particle 

F;{x, k ± ;x', k' ± ) := F 3 (l - x, -kr, 1 - x' , -k' x ). 



If the function additionally depends on the component of 
the total angular momentum J z = n, a tilde symbolizes 
the operation Fi{n) = n). 



1. Helicity table of the annihilation graph 



s z = s'= 0. 



(A5) 



This means that the seagull graph does not contribute 
when J z ^0, because it has a factor proportional to 
cos(</3 — ip') resulting from (A5). A rather surprising con- 
sequence is that the dynamic graph contributes only for 
J z = ±1, the instantaneous graph only for J z — 0. How- 
ever, to maintain rotational invariance, both diagrams 
must yield the same value, though one shows much more 
structure than the other. Degeneracy of the orthopositro- 
nium ground state with respect to J z is found without 
the annihilation channel and inclusion must not destroy 
it. At hrst glance, there seems to be a manifest breaking 
of rotational symmetry: the helicity table |v| separates 
between states with (Ai,A2) = (TT) and (Ai, A 2 ) = (||). 
But this is only a consequence of the integration over 
the angles: for J z = +1 the (J,!)- combination gives no 
contribution, and likewise does the (1"f)-combination for 
Jz - -1. 



The functions Fi{X;2) :— Fi(x, k±; x' , k' ± ) displayed in 
Table are: 



Fi(x, k±:x', k'j_) 



a 2 



J* l.l 



TT LU* XX'(1 — X)(l — X 1 ) 

F 2 (x, k ± ;x', k' ± ) := ^ ( + 43 \ 

TT \U> XX ) 



F 3 (x,k±;x',k' 1 _) 
F 4 (x,k ± ;x',k' ± ) 



a 2 m k'^ 

* Ai— jz 7<5|^ |,i 

TT UJ* X(l — X) 1 — X 

2 fcxfcj 

TT \U) X (1 — X) 



The table for J z — — 1 is obtained by inverting all he- 
licities. Note that the table has non-vanishing matrix 
elements for \J Z \ < 1 only. This restriction is due to the 
angular momentum of the photon. 
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out: in 


AiA 2 =n 


AiA 2 =U 


AiA 2 =1T 


(Ai,A 2 =|| 


AiA 2 =TT 


^i(i;2) 


F 3 (2;l) 


^3*(2;1) 





AiA 2 =n 


F 3 (l;2) 


F 2 *(l;2) 


^4(2;1) 





AiA 2 =4T 


F£(i; 2) 


fk(l;2) 


F 2 (l;2) 





AiA 2 = Lj 















TABLE V. Helicity table of the annihilation graph for 
J z > 0. 
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FIG. 10. Comparison of multiplets for n=2: (a) results of the present work with a — 0.3, A = 1.0 m/, Ni = N 2 = 21; (b) 
equal-time perturbation theory (ETPT) up to order 0(a 4 ). 
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